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ABSTRACT 

It has been widely thought that measuring the misalignment angle between the 
orbital plane of a transiting exoplanet and the spin of its host star was a good discrim- 
inator between different migration processes for hot-Jupiters. Specifically, well-aligned 
hot- Jupiter systems (as measured by the Rossitcr-McLaughlin effect) were thought to 
have formed via migration through interaction with a viscous disk, while misaligned 
systems were thought to have undergone a more violent dynamical history. These con- 
clusions were based on the assumption that the planet-forming disk was well-aligned 
with the host star. Recent work by a number of authors has challenged this assumption 
by proposing mechanisms that act to drive the star-disk interaction out of alignment 
during the pre-main sequence phase. We have estimated the stellar rotation axis of a 
sample of stars which host spatially resolved debris disks. Comparison of our derived 
stellar rotation axis inclination angles with the geometrically measured debris-disk 
inclinations shows no evidence for a misalignment between the two. 
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1 INTRODUCTION 

The discovery of planets beyond the confines of our Solar 
system has presented many surprises and continues to chal- 
lenge our understanding of planet formation and their subse- 
quent evolution. This is particularly true in the case of hot- 
Jupiters, whose short orbital periods of a few days or less 
was unexpected - under the standard core-accretion theory 
of planet for mation, volatile gas- giants should form beyond 
the snow- line (|Pollack et al.|[l996l h It is now widely accepted 
that hot-Jupiters did not form in-situ at their current loca- 
tions, but that some mechanism caused their inwards mi- 
gration towards their parent star. 

A number of theories have been postulated to ex- 
plain planetary migration. One possible mechanism for 
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forming short-period gas-giants is the pumping of ini- 
tially wide circular orbits to high ecce ntricities. This could 
occur via planet-planet scatte ring (|Rasio k, Fordl 1 19961 : 
IWeide nschilling & Marzaril Il996l ), or perturbations from a 
distant stellar binary companion (|Eggenberger et al]|2004l '). 
The highly eccentric orbit then brings the gas-giant suffi- 
ciently close to the host star that tidal dissipation quickly 
draws the planet to a new, smaller orbital separation. In 
this scenario, the interactions and scattering involved may 
lead to large changes in the value of the orbital inclination. 
Interactions between the planet and a viscous disk, on the 
other-hand, may also drive the planet inwards but is not 
thought to perturb the initial orbital inclination. 

The close alignment of the rotation and orbital axes 
in the Solar system (~7°; lBeck fc Giledl2005l ) is attributed 
to the formation of the Sun and planets from a single ro- 
tating proto-stellar disk which was also initially coplanar 
to the solar-rotation axis. On the premise that disks and 
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stellar rotation axes are aligned, Rossiter - McLaughlin (RM) 
observations of transiting systems (e.g. iTriaud et al.ll201ol 
and references therein) have sought to discriminate between 
migration caused by planet-disk interactions (leading pre- 
sumably to aligned systems) , and migrations involving some 
violent dynamical history (leading to misaligned systems). 
Recent theoretical work, however, challenges the view that 
the stellar spin axis and the disk rotation axis should be 
aligned. For example, in numerical sim ulations of star for- 
mation, iBate. Lodato Sz Pringld (|2010h found that the rota- 
tion axis of the final disk may be heavily governed by the 
angular momentum component of the material that was last 
accreted. While this can lead to a misalignment between the 
disk and star, it is only thought to be significant for light 
disks which may have in sufficient mass to form giant planets. 

On the other-hand. lLai et al.l (|2010l ) present arguments 
that the observed star-orbit misalignment could instead re- 
sult from alterations in the stellar spin axis, introduced by 
the star- disk interaction durin g the pre-main-sequence phase 
(also see iFoucart fc LaH boiO,). Lai et al. consider the well 
known fact that a magnetic protostar exerts a warping force 
on th e inner part of the accretion disk (e.g. iBouvier et al.l 
l2007h . Previous authors have assumed that this results in 
significant warps to the inner disk, whereas lLai et al] (|2010l ) 
argue that viscous processes in the disk itself will smooth 
these torques, resulting in a largely unwarped inner disk. 
Given a flat disk, the torques arising from the star-disk in- 
teraction will act on the star itself, changing the stellar spin 
axis on a timescale given by 
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where M* and R* are the mass and radius (in solar units) of 
the protostar, respectively, M is the accretion rate in solar 
masses per year, n n is the inner radius of the accretion disk 
in stellar radii, u s is the spin rate of the protostar and Q(ri n ) 
is the rotation rate of the accretion disk at the inner disk 
radius. 

However, the mechanism proposed by lLai et all (|20ld ) 
may not be effective in practice, as the timescale for spin 
evolution, t 3 pi n , is of the same order as the disk evolu- 
tion timescale. Near- infrared observations of protostars show 
that the majority of protostellar disks h ave dispersed by 
the age of 5 Myr (|Hernandez et all I2008T ) whilst observa- 
tions of the accretion rates onto young stars also show that 
the accretion rate declines r apidly with increasing age and 
decreasing stellar mass fe.g ISicilia-Aguilar et al.ll2006l ). so 
the accretion rate on many protostars may well b e below 
the canonical 10 _s M^yr -1 assumed by Lai et al.l (|2010l) - 
In addition, the results of lLai et al.l l|201(j ) rely on the in- 
ner disk not being 'significantly warped', however there is 
good evidence that the inner disks of some young stars do 
contain significant di sk warps (see e.g. IBouvier et al.ll2007l ; 
iMuzerolle et alJl200St ). 

If correct, these theoretical papers potentially have im- 
portant ramifications for our interpretation of the results of 
RM observations. Indeed, if the stellar rotation axis can be 
driven from coplanarity with the surrounding disk, or vice- 
versa, then RM observations would essentially be rendered 



useless as a tool for determining the migration mechanism 
responsible for forming hot- Jupiter's. For these reasons, it is 
important to seek observational evidence for such processes. 
In this paper we present a study of star-disk alignment in 
debris disk systems. 



2 MEASURING THE STAR-DISK 
ALIGNMENT 

For the purposes of this work, we have concentrated on sys- 
tems with spatially resolved debris disks. We have also as- 
sumed that the debris-disk plane is representative of the pri- 
mordial disk and, likewise, that the presently observed stel- 
lar orientation is the same as the protostars'. The inclination 
of the disk to our line-of-sight can then be measured geomet- 
rically by calculating the fore-shortening of the semi-minor 
axis of the disk relative to the semi-major axis (although in 
reality the models used to determine the disk geometry are 
somewhat more complex). 

A more indirect approach is needed in order to de- 
termine the inclination angle of the stellar rotation axis, 
however. To do th is we have followed the method of 
IWatson et all (|2010l ) who compiled the stellar rotation in- 
clination angles for 117 exoplanet host stars, and we refer 
the reader to that paper for in-depth details of the meth- 
ods used, as well as a discussion on possible sources of sys- 
tematic errors inherent in the technique. In summary, it is 
possible to determine the inclination angle, i, between the 
rotation axis of a star and the observers line-of-sight from 
measurements of the projected equatorial velocity (usini), 
the stellar rotation period (P ro t) and the stellar radius (R*) 
via the equation 
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The projected equatorial rotation velocity, t>sini, can be 
measured using high-resolution spectroscopy, while the stel- 
lar radius can also be indirectly determined from spectra 
or, less frequently, dir ectly via interferometr y, lunar occul- 
tations or eclipses (e.g. lFracassini et al.ll200iT ). Precisions on 
stellar radiu s measurements of ^3 pe r cent are now regularly 
quoted (e.g. iFischer fc ValentilboOSl 1 ). 

Determining the stellar rotation period, on the other- 
hand, tends to be more troublesome. For some active stars, 
the stellar spin period can be determined photometrically to 
high precision by tracking the passage of large star spots on 
their surfaces. For those systems which do not have photo- 
metrically measured rotation periods, measurements of Ca 
n H and K emission can be used to estimate the rotation pe- 
riod by applying t he chromospheric e mission - rotation pe- 
riod relationship of lNoves et al.l (| 19841 ) . Naturally, this latter 
method is less precise, and is also affected by intrinsic vari- 
ability of the Ca II H and K emission due to, for example, 
solar-like activity cycles or the rotation of magnetic regions. 

We have carried out an extensive literature search and 
present vsmi, R*, and P ro t estimates for a number of main- 
sequence stars which host spatially resolved debris disks in 
Table [T] Since one of the pre-requisites for measuring a 
stellar rotation period is that the star must be magneti- 
cally active, we are restricted to lower main-sequence stars 
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later than ~F5V which have a convective envelope (and 
are thereby capable of sustaining a stellar dynamo). Of the 
20 main-sequence stars with resolved debris disks, only 10 
have spectral types of F5V or later. Of these, we can find 
no recorded Ca II H and K emission measurement for HD 
181327, and is therefore omitted from our list. 

We should note that we have not considered pre-main 
sequence stars in our analysis. This is for two principal 
reasons. First, given their fully convective nature, it is 
not certain t hat th e activity-rotation period relationship of 
iNoves et all (|l984T l (which was calibrated for main-sequence 
stars) holds, indeed an entirely different stellar dynamo 
mechanism may o perate in pre-main sequence stars (e.g. 
IScholz et al]|2007f ). Second, radius estimates for pre-main 
sequence stars are also notoriously unreliable, since they de- 
pend upon age estimates which are uncertain by a factor 
of se veral (e.g. iNavlor] [2009; Bar affe. Chabrier fc Gallardol 
120091 ). 

2.1 Adopted stellar parameters and errors 

In order to determine sini via equation [2] we have taken a 
weighted mean of the entries in Table [1] for the final values 
of vsini and R, . In id entical fashion to that carried out 
in IWatson et ail j201Ch . where no error was quoted on a 
published vsini value we have taken it to be 1.0 kms -1 
(twice the typical er ror assumed on vsini m easurements, 
see the catalogue of iFischer fc Valentil 120051 for example). 
Regarding published radii with no associated error bar, we 
have taken the error to be 10 or 20 per cent of the absolute 
value. The choice between 10 or 20 per cent is taken to 
ensure that radii estimates with associated error bars were 
given a higher weighting than those without formal errors. 

For stars with photometrically derived rotation periods 
which have no associated error bar, we have taken the error 
to be 10 per cent. This is commensurate with the typical 
error bars quoted on such measurements. Where available, 
photometrically derived rotation periods are adopted, oth- 
erwise the rotation period is estimated from th e strength of 
the C a H H and K emission (|Noves et alJll984h . Again, fol- 
lowing IWatson et all l|2010h . for each log R' HK measurement 
reported in Table[T]we have determined, where possible, the 
number of observations and period over which they were 
carried out (see Table HJ. Where details are not present, 
or are ambiguous, we have assumed they are from a single 
observation a n d hav e flagged them as 'individual?'. As in 
IWatson et all (2010), each star was assigned a grade of P 
(Poor), O (O.K.), G (Good) or E (Excellent) based on how 
well monitored it was. We then assigned general error bars 
on the logi?^^ values dependent on their assigned grades 
and spectral type. These error bars are derived from the av- 
erage rotationally modulat ed variations outlined in Section 

3.1 of lWatson etafl (|2010l ). For a detailed discussion of the 
systematic errors on the derived parameters, we refer the 
reader to this work. Table [2] lists the adopted parameters 
for each star. 

2.2 Determining the stellar inclination angle 

Equation [2] can be thought of as a naive estimator of sin i 
as it is geometrically unconstrained (e.g. sini > 1 is al- 
lowed). While a value of sini > 1 is unphysical, it does 



allow potential prob l em ca ses to be identified. Again, we 
follow Watson et al l j2010h and reject systems with sini's 
that are 1-a greater than 1 from further analysis - flagging 
these as having a high probability of being affected by sys- 
tematic errors. This results in the omission of 2 systems, 
HD 53143 and HD 139664, both of which have naive sini 
estimates significantly greater than 1 (see the first two en- 
tries of Table [2]). In the case of HD 139664, the B - V value 
places it at the extreme edge of th e chromospheric e mis- 
sion - rotation period calibration bv lNoves et alj (| 1984T ) . In 
addition, the star is classified as having a luminosity class 
IV, and therefore both the derived rotation period from the 
INoves et alj (1 19841 ') relationship (which is only calibrated for 
main-sequence stars) and radius may also be suspect. HD 
53143, on the other hand, is more problematic. It appears 
to have a secure rotation period which has been measured 
photometrically and that also agrees very well with the pe- 
riod derived from the Ca II H and K emission. In addition, 
all of the measured v sini's and radii are consistent with one 
another. Yet, despite this and the fact that it appears to be 
a solid main- sequ ence star with an age of 1.0 ± 0.2 Gyr 
(|Kalas et alj|200a ). we derive sini ~ 1.5 ± 0.4. We can only 
assume that 1 or more of the measurements are affected by 
systematics. 

For the 8 remaining systems we have carried out a 
Markov-chain Monte Carlo (MCMC) analysis which not 
only provides a means of optimising the fit of a model to 
data but explores the joint posterior probability distribution 
of the fitted parameters and allows proper 1-a two-tailed 
confidence limits to be placed on the derived sini's. In addi- 
tion, MCMC rejects unphysical combinations of parameters 
that result in sini > 1. For the purposes of this work, we 
have f ollowed the MCMC process outlined in I Watson et al.l 
(|201Ch . keeping the same 1000-step burn-in phase and carry- 
ing out 1,000,000 jumps. The results of this MCMC analysis 
are shown in Table [5] 



3 RESULTS AND DISCUSSION 

Table [3] shows our derived stellar rotation inclination angles 
versus published debris-disk inclinations. As can be seen, 
there is no obvious evidence for large mis-alignments of the 
stellar rotation axes and debris-disk planes in any of these 
systems. By the nature of the method, the best constrained 
systems have sini ~ 0.5 (i, =30°). This is because at high 
inclinations the sine curve is relatively flat, and thus small 
errors in sini (which is what is directly calculated from the 
observables in equation [2| propagate to form large errors 
when expressed in degrees. At sini's of ~0.5, the sine curve 
is much steeper, and travelling along the sine curve does 
not vary the inclination i, as quickly as it does at high 
sini's. As one moves to lower sini's, measurement errors on 
vsini naturally increase as the projected rotational broad- 
ening decreases. The fact that the best constrained systems, 
HD 22049 and HD 107146, with errors on i» of only 5-9° 
appear to align closely with their debris disk gives us both 
confidence in the technique, and further strengthens our as- 
sertion that we see no evidence for a detectable difference 
between the sky-projected angle of the disk and the that of 
the stellar rotation axis. In addition, it should be noted that 
HD 22049 is known to host a planet that has had the inclina- 
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tion of its orbital plane accurately determined to be ipianet = 
30.°1±3.°8 - s uggesting coplanarity between the planetary 
orbit and disk ijBenedict et alJbood ). Furtherm ore, star spot 
model ing of a MOST light curve of HD 22049 bv lBrvce et ail 
(2006) determined the inclination of the stellar rotation axis 
to be i* = 30° ±3°, in excellent agreement with our derived 
values. We do caution, however, that the absolute direction 
of the axis (whether the rotation axis is pointing towards or 
away from the observer) cannot be ascertained, and there- 
fore we do not have a knowledge of the full three-dimensional 
geometry of the star-disk systems. 

We can test the significance of our result using a rank- 
order approach. If the disk and stellar spin axes are closely 
aligned, we expect a ranking by disk inclination to agree 
well with a ranking by stellar inclination. This is indeed 
what we find. We calculate the spearman rank-order corre- 
lation coefficient for our data, and find a value of 0.82. If 
we repeat this analysis for all 40320 random permutations 
of our stellar spin data, only in 0.1% of cases do we find a 
better correlation between stellar spin and disk inclinations. 
This suggests we reject the null hypothesis (that stellar spin 
and disk inclinations are uncorrelated) at a significance of 
99.9%. 

We performed a further test of significance by first as- 
suming that the stellar spin and disk axes were uniformly 
distributed around the sky. We drew 8 values of disk and 
stellar inclinations from a uniform distribution in cosi, and 
took account of our errors by perturbing these inclinations 
by an amount equal to the errors on our observations. For 
the inclinations, which have asymmetric errors, we used an 
average of the two errors. For the disks we assigned an er- 
ror of 3 degrees where none was quoted. Where the given 
disk inclination is a range, we set the error to half of that 
range. We computed 1 x 10 6 such simulated datasets, and 
only 0.4% showed a better correlation (as judged by the 
Spearman rank-order coefficient) than our original data. We 
therefore reject the hypothesis that the stellar spin axes are 
independent of the disk inclination axes with a significance 
of 99.6%. 

In reality, it is likely that the stellar spin axes and disk 
axes are aligned to some degree of precision. It is reason- 
able to ask what is the largest degree of misalignment per- 
mitted by our data. Answering such a question would re- 
quire a full Bayesian treatment which accounts for the fact 
that we only observe the component of misalignment along 
the line of sight. Such an analysis is quite complex (see 
iFabrvckv fc Wirmll2009l . for example), and is left for a future 
work. 

A recent an a lysis of Rossiter-McLaughlin observations 
bv iTriaud et all lE)lch suggest that between 45 - 85 per 
cent of hot-Jupiters appear to be significantly misaligned. 
Our work in this paper reveals no similar degree of misalign- 
ment between debris disks and their host stars. We note, 
however, that all of the systems in our study have host stars 
with e ffective temperature s below ~6140K (see Table[3]). Re- 
cently, IWinn et all (|2010h have highlighted that exoplanet 
host stars with effective temperatures below ~6250K ap- 
pear to have the planet-star spin axes preferentially aligned, 
whereas exoplanets orbiting hotter host stars are more likely 
to be misaligned. It is, therefore, possible that a yet to 
be determined mechanism which only drives star-disk mis- 
alignments in hotter systems could be operating which our 



small study has missed. We conclude, however, that there 
appears to be no substantial evidence to s uggest tha t a uni - 
vers al process, such as that outlined by lLai et all (l20ld ) 
and lBate et"afl l|201Cl ). is a major mechanism in misaligning 
planetary orbits. 
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Table 1: Published data on the properties of 10 stars hosting resolved 
debris disks. Rotation periods quoted with no refer ence have been cal - 
culated using the adjacent \ogR' HK entry using the iNoves et all l| 19841) 
chromospheric emission - rotation period relationship along with (B — V) 
values taken from NStED. 
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Table 2: Adopted parameters and sini estimates for the stars in our 
study. The colour index (B — V) is only indicated for stars where the 
rotation period has been deter mined from t he ch romospheric activity - 
rotation period relationship of lNoves et al, (1984). The first two entries 
have 'naive' sini estimates in column 9 (indicated by an asterisk) as 
derived from equation [2] which, complete with their formally propagated 
errors, result in sini's significantly above 1. These stars are omitted from 
further analysis. For the remaining eight stars, column 9 gives the final 
derived sini value, followed by the l-o two-tailed confidence limits, as 
derived from a Markov-chain Monte Carlo analysis. 
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8.599 


1.000 


5.419 


2.108 


0.829 


0.048 


0.742 


0.999 


0.123 


0.000 


92945 


5.022 


0.468 


7.176 


2.830 


0.786 


0.024 


0.894 


0.908 


0.091 


0.087 


107146 


5.000 


0.447 


3.496 


1.35 


0.993 


0.014 


0.611 


0.353 


0.141 


0.138 


197481 


8.832 


0.959 


4.846 


0.20 


0.835 


0.018 




0.999 


0.062 


0.000 


207129 


2.319 


0.447 


17.129 


1.610 


1.048 


0.018 


0.600 


0.746 


0.167 


0.187 
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Table 3. Comparison of the derived stellar rotational axes and 
published disk-plane inclinations. For HD 10647 and HD 10700 
the lower value for the disk inclination corresponds to that derived 
from the observed disk dimensions and which we take to be the 
most probable value. References for the disk inclinations are given 
in the fourth column. Estimates of the host star mass and effective 
temperature from the NStED database are quoted in the final two 
columns. 



HD u (°) idisk (°) ref. Mass (M Q ) T eff {K) 



10647 


49111 


^52 


fLiseau et al. 2008) 


1.20 


6140 


10700 


4&±% 


60-90 


f Greaves et al. 2004) 


0.91 


5500 


22049 


31±I 


25 


f Greaves et al. 1998) 


0.78 


5090 


61005 




SO 


(Maness et al. 2009) 


0.89 


5440 


92945 


65+ 21 
D0 -10 


70 


fKrist et al. 2005) 


0.77 


5060 


107146 


2lig 


25±5 


fArdila et al. 2004) 


1.09 


5850 


197481 


90l» 


90 


CKrist et al. 2005) 


0.49 


3560 


207129 


47t a J 


60±3 


fKrist et al. 2010) 


1.11 


5890 
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Table 4: Compilation of chromospheric indices (\ogR' HK ) for the stars in 
Table 1 for which no measured rotation periods have been reported. The 
spectral type of the host star is given in column 2. Entries in bold give the 
grade assigned to each star (P = Poor, O = O.K., G = Good, and E = 
Excellent) followed by the weighted mean of the log R'hk measurements 
and adopted error bar (see section r2.ll for details). Reference numbers 
are identical to those used in Table [1] 



Name 


Type 


log K HK 


Observations 




Kef. 


ttt\ i r\r< A 

HD 10647 


F8V 


-4.680 


individual? 




8 






-4.700 


individual on 2001 Aug 04 




9. 






-4.714 


individual? 




11. 








(P) Adopted value: -4.698 


± 0.060 




HD 61005 


G3/5V 


-4.260 


1 obs on UT 14/12/1992 




5 






-4.324 


individual? 




11 






-4.360 


1 obs on 28/10/2002 




15 






-4.337 


individual? 




16 








(P) Adopted value: -4.320 


± 0.075 




HD 92945 


K1V 


-4.320 


13 obs in 6 months. Report a 


= 2.72% 


3 






-4.393 


individual? 




16 








(O) Adopted value: -4.325 


± 0.077 




HD 107146 


G5 


-4.340 


8 obs in 5 months. Report a = 


= 3.04% 


3 








(O) Adopted value: -4.340 


± 0.057 




HD 139664 


F3/5V 


-4.621 


individual? 




11 








(P) Adopted value: -4.621 


± 0.060 




HD 207129 


GOV 


-4.800 


1 obs on UT 28/06/1993 




5 






-4.850 


1 obs on 2004 Aug 23/24 




9 






-5.020 


individual? 




16 



(P) Adopted value: -4.89 ± 0.075 
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